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ABSTRACT 

We use stellar dynamics, strong lensing, stellar population synthesis models, and weak lensing shear 
measurements to constrain the dark matter (DM) profile and stellar mass in a sample of 53 massive 
early-type galaxies. We explore three DM halo models (unperturbed Navarro Frenk & White [NFW] 
halos and the adiabatic contraction models of Blumenthal and Gncdin) and impose a model for the 
relationship between the stellar and virial mass (i.e., a relationship for the star-formation efficiency 
as a function of halo mass). We show that, given our model assumptions, the data clearly prefer a 
Salpeter-likc initial mass function (IMF) over a lighter IMF (e.g., Chabrier or Kroupa), irrespective 
of the choice of DM halo. In addition, we find that the data prefer at most a moderate amount of 
adiabatic contraction (Blumenthal adiabatic contraction is strongly disfavored) and are only consistent 
with no adiabatic contraction (i.e., a NFW halo) if a mass-dependent IMF is assumed, in the sense of 
a more massive normalization of the IMF for more massive halos. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: halos — galaxies: structure — 
dark matter 



1. INTRODUCTION 

Numerical simulations of cold dark matter (CDM) ha- 
los predict density distributions with an inner logarith- 
mic slope dlnp/dlnr = — 1 and a slope o f dlnp/dhar = —3 
at large radii (|Navarro et aljl99ftl20lcl hereafter NFW). 
However, the inclusion of baryons - and the dissipa- 
tional processes associated with them - in galaxy for- 
mation simulations may significantly alter the central 
dark matter distribution. The DM may show an in- 
creased central density due to radiative dissipation from 
the baryons that depends on the mass of stars formed in 
the g alaxy (e.g., the adiabatic contraction, or AC , mod- 
els of Bl umenthal et all 1198ft iGnedin et all 12004 here- 
after B86 and G04) or dynamical friction may counter- 
balance this effect (Abadi et al. 200 9]) and could result 
in cored profiles (e.g.. lEl-Zant et all 12004 iNipoti et all 
12004 ). The form of the CDM halo is critical for under- 
standing the relative importance of baryonic and total 
mass in governing the structure of early-type galaxies 
(ETGs). 

Strong lensing can measure masses to a few percent 
precision and is therefore a particularly powerful probe 
of this regi me, by itsel f and in combination with other 
techniques (jTreul l2010l and references therein). Stud- 
ies have previously used strong and weak gravitational 
lensing (SL and WL) to constrain the distribution of 
DM for massive ETGs and found that the data are 
broadly and on average consistent with a NFW halo (e.g., 
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IGavazzi et all 120071 ) . However, the SL and WL analy- 
sis alone cannot simultaneously probe the stellar mass 
M„ and central dark matter slope due to degeneracies 
(|Lagattuta et al.ll2009D : additional information must be 
included to disentangle the contributions of stars and 
CDM to the central density distribution (s imilar to the 
disk-h alo degeneracy for spiral galaxies). ISchulz et all 
(2009) have recently used WL and stellar velocity dis- 
persions to investigate evidence for AC in SDSS ETGs. 
They find that the stellar velocity dispersions require 
more central mass than would be inferred from the stellar 
mass (assuming a Kroupa initial mass function, or IMF) 
and the central DM mass extrapolated from a NFW halo 
fit to the WL data. However, their findings are also con- 
sistent with no AC if a Salpeter IMF is assumed for the 
st ellar component. 

Jiang & Kochanek (2007) point out that the degen- 
eracies between the central CDM slope and M* can 
only be broken by using at least three mass probes. 
They use SL and dynamics to probe the central bary- 
onic and CDM distribution of ETGs and include an en- 
semble measurement of the halo mass from WL data 
to find th at AC models are f avored over a NFW halo. 
However, Uiang fe' Kochanck (2007J) assume a redshift- 
dependent but otherwise constant stellar mass-to-light 
ratio M*/L for each of the galaxies, althoug h the IMF 
(and therefore M*/L) may be non- universal (jTreu et all 
2010). A full understanding of the central mass distribu- 
tion of galaxies requires constraints on both the stellar 
and CDM components as a function of mass. Although 
this joint analysis is more complex and there are resid- 
ual degeneracies, it provides a rare opportunity to con- 
strain directly the IMF of distant massive ETGs where 
resolved stellar population diagnostics are not avai lable 
(see, e.g.. lCappellari et alll2009l : iBastian et al.ll2010L and 
references therein). 

In this Letter we present simultaneous constraints on 
M* and the form of the central dark matter distribution 
for a sample of massive ETGs. We use SL, stellar dy- 
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namics, stellar population synthesis (SPS) stellar mass 
estimates Mf FS given an IMF, and WL shear measure- 
ments to constrain a model for the relationship between 
stellar and virial mass. We allow for several different DM 
halo profiles and we fi t for an IMF mis match parameter 
relating M* to M^ PS (jTreu et al.ll2010D. The two m ain 
improvements over the analysis bv lTreu et al.l (J2010D are 
the inclusion of WL data and the adoption of more flex- 
ible models. We use a standard flat ACDM cosmology 
with Hm = 0.3 and h = 0.7. 

2. THE DATA AND MODEL 

Our sample is derived from the Sloan Lens ACS Sur- 
vey (SLACS: IBolton et al.l [20061: IGavazzi et al.1 [20071: 
IBolton et aLll2008H Auger et al.ll2009H . which consists of 
nearly 100 strong gravitational lensing galaxies. We ex- 
clude 6 galaxies which appear to be structurally dif- 
ferent than the re st of the SLACS early-type lenses 
(|Auger et alJ^OlOal ). These systems are significant out- 
liers from the otherwise very tight relation between ef- 
fective radius, velocity dispersion, total central mass, 
and stellar mass; this is likely due to systematic errors 
in t he stellar velocity disper sions of these systems (also 
see IJiang fc KochaneH I2007D . We therefore use a sub- 
sample of 53 ETGs that have well-determined central 
stellar ve l ocity dispe rsions, SPS stell ar masses assuming 
iChabrierl ([200l and iSalpeterl (fl955l) IMFs, and central 
projected mass estimates from SL. We also use WL shear 
measurements for a subset of 44 lenses with deep ACS 
imaging. The details of how th ese data a r e derived can 
be found elsewhere (IG avazzi et al. 2007; B olton et al.l 
120081 lAuger et all 120091 l2010af K We emphasize that 
the SLACS lens galaxies are indistinguishable from twin 
ETGs selected from SDSS to have the same velocity dis- 
persion and redshift and can therefore be considered as 
representative of the overall population of massiv e ETGs 
(ITreu et al.N2006L [2001 [Auger et al.H200l I2010afl . 

We model the lenses as two-component mass distribu- 
tions consisting of stars and a CDM halo. Both com- 
ponents are modeled as spherical since our data do not 
constrain the halo flattening and our dynamical analy- 
sis is restricted to spherical Jeans modeling. Compar- 
isons with non-spherical models shows that the most 
relevant quantities such as mass-density slope and to- 
tal m ass are relatively ins ensitive to this approximation 
(e.g.. lBarnabe et al.|[2009h . 

The stars are modeled with the IHernquistl (1990) pro- 
file and the scale radius is set by the rest-frame y-band 
effective radius, r a — 0.551^; our conclusions are un- 
changed if we adopt a Uaffd (|1983D profile. We assume 
isotropic orbits, consistent with typical findings for mas- 
sive E TGs (Barnabc et al. 2009). Mild radial anisotropy 
(|Koopmans et al.1 12009') only alters the inferred velocity 
dispersions by a few percent, i.e. an amount smaller than 
the errors, and therefore does not change our conclusions. 

We investigate three different density profiles for the 
CDM halos. First, we consider the NFW profile with 
the concentration fixed by c(M v ; r ) oc M"^ 094 follow- 
ing IMaccio et al.l ((2P08i) , neglecting intrinsic scatter for 
computational simplicity. We also modify the NFW 
halo to include AC using the prescriptions of B86 and 
G04, with the baryon fraction set by Equation 1 and 
the scale radius given by r a ; these three profiles are 
the most widely adopted descriptions for CDM ha- 



los (e.g.. IJiang fc Kochanekl [20071 iSchulz et al.l [200l 
iNapolitano et al.H2010D . and the NFW (B86) halo has 
the lowest (highest) central DM density. 

We fix the relationship betw een stellar an d viria l mass 
to have the form suggested bv lMoster et al.l (|2010f) (who 
use a Kroupa IMF) based upon abundance matching 
techniques. We assume 



M» =Mf PS =M* - 



(M vir /M viri0 



,7i 



1 + (M vir /M vir , r 



(71-72)/^ 



(1) 



where Mf PS is the stellar mass inferred from the SPS 
models assuming a Salpeter IMF, 71 = 7.17, 72 = 0.201, 
and P = 0.557 and we fit for M* j0 and M v j ri0 . We note 
that our results are qualitatively the same if we use dif- 
feren t values for 71,72, and (3 (e.g., from IDutton et al] 
l2010f ). although the goodness-of-fit changes; we will ex- 
plicitly fit for these parameters in a future paper to as- 
certain the extent to which one model prefers a given set 
of parameters over another. 

We also allow for an offset between the model stellar 
mass and the SPS stellar mass, equivalent to a mass- 
dependent IMF (referred to as the 'Free' IMF model). 
The offset has the form 



log 



Mjps 
M* 



= -??log 



M. 



10 n M^ 



log a, 



(2) 



where a is equivalent to that introduced by ITreu et al.l 
(|2010() and rj is equivalent to the slope of Equation 
5 i n ITreu et al.l (I2010D . We adopt the Mf s derived 
by I Auger et al.l (|2009f ). although our results are robust 
with respe ct to changes in the SPS models and choice 
of priors (jTreu et al.l 12010ft due to the relative sim- 
plicity of the older stellar populations of the SLACS 
ETGs. Uncertainties in SPS models of evolv ed galaxies 
(|Conrovfc Gunnl 120101 : iMaraston et~aTl 120091 ) may lead 
to systematic off sets in Mf PS of order ~ 0.05 dex (e.g., 
ITreu et al.l[20Tol) . 

We perform a Markov Chain Monte Carlo simula- 
tion to fit the models to the data; we propose M v ; r for 
each lens and use Equation 1 (and 2 for the Free IMF 
model) to predict a stellar mass. These masses normal- 
ize the halo and bulge, and we use these profiles to pre- 
dict the stellar kinematics, SL, and WL signals for each 
galaxy. We compare these quantities and the model stel- 
lar masses with the observed quantities on a galaxy-by- 
galaxy basis (i.e., we do not bin or average measurements, 
including the WL) to optimize the free model parameters 
and determine goodness-of-fit values. 

3. RESULTS 

The best-fit parameters and x 2 for each model are 
listed in Table [T| and we show the best- fit mass distribu- 
tions for a characteristic lens system (Ly = 1O 1O85 L0) 
in Figure [T] This illustrates that the low normalization 
of the Chabrier IMF requires a more massive CDM halo 
to increase the central projected mass and be consistent 
with the mass required by the SL data (this is true for 
the NFW and AC models). However, increasing the halo 
mass also increases the scale radius because the virial 
radiu s increases and the concentration becomes smaller 
(e.g.. IMaccio et aLll2008t ): therefore, the halo mass grows 
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at a faster rate than the projected dark matter mass in- 
side the Einstein radius. The net effect is that very mas- 
sive halos are needed to match the SL data if a Chabrier 
IMF is used (see Figure [2]). 

These massive halos cannot fit the SL masses and the 
WL shears (the inferred M* values are also much higher 
than the measurements), and the data therefore clearly 
disfavor a Chabrier IMF (Table[T]), consistent with the re- 
sults from SL and dynamical constraints only (jTreu et al.l 
120101 ). The G04 model is slightly preferred over the NFW 
halo for the Salpeter IMF while the B86 model is disfa- 
vored. The normalization and break point of the M v ; r -M* 
relation are consistent between the NFW and G04 model 
with a Salpeter IMF, and bo t h are also consistent with 
the results from lMoster et all (|2010f ) when the ~ 0.2 dex 
offset between Kroupa and Salpeter IMFs is taken into 
account. 

Figure [2] illustrates that the favored models require sig- 
nificantly higher star-formation efficiencies compared to 
the values of « 0.02 - 0.1 (for a Chabrier IMF; double 
this for Salpeter) found for massive ETGs from galaxy- 
galaxy WL studies (e.g., Mandclbaum et al. 2006), which 
may point to a generic inadequacy of the family of mod- 
els considered here. Interestingly, simulations of galaxy 
formation may require similarly high star-formation ef- 
ficiencies to reproduce the observed isothermal centra l 
mass distributions of massive ETGs (|Duffv et al.ll2010l) 
and suggest that the high efficiency may not be due to 
our simplifying assumptions. 

The results for the free IMF model are somewhat more 
difficult to interpret due to the addition of two extra pa- 
rameters; the fit is significantly improved compared to 
the fixed IMF models, as is expected when allowing for 
a more flexible model. The NFW model and G04 model 
fit the data equally well, but again the B86 model is dis- 
favored. Note, however, that the form of the M V j r -M* 
relation that we are using was derived assuming a con- 
stant IMF; allowing for a varying IMF would necessarily 
change the structural parameters 71,72, and /?. Never- 
theless, the present model indicates that even the AC 
halos require the IMF to have a slight trend with mass 
in the sense that more massive galaxies require more 
Salpeter-like IMFs. The trend is stronger with the NFW 
ha lo (Figure B t and is consistent with the relation found 
bv lTreuet airpoToll . 

4. DISCUSSION 

There is a growing consensus that the fraction of 
DM at the center of massive ETGs increases with 
galaxy mass dTortora et al.ll2009b INapolitano et al.ll20Tot 
iGraves fe Faberll2010HAuger et alJl2010aj ). This increase 
could be due to a genuine increase in cold dark mat- 
ter, perhaps as the result of AC me diated by the de- 
tailed accretion history o f halos (e.g., lAbadi et al.l 120091 : 
lLackner fc Ostrikerll2010H . Alternatively, the trends may 
be the result of a mass-dependent IMF leading to more 
baryonic DM (in the form of low-mass stars or stellar 
remnants, depending on whether the IMF is bottom- 
heavy or top-heavy, respe ctively) in more massive galax- 
ies fe.g.. lTreu et al.ll2010f). perhaps due to cosmic evolu- 
tion of the IMF (|van Dokkum et all 120081 ) and a mass- 
dependent formation redshift for the stellar mass (e.g., 
iThomas et aDEQQJ Ivan der Wei et al.ll2009h . 



4.1. Ruling Out "Light" IMFs for Massive ETGs 

We find that our sample of massive ETGs is inconsis- 
tent with a Chabrier IMF, even after allowing for strong 
CDM halo contraction and a mass-dependent IMF. In- 
deed, our data disfavor common "li ght" IMFs gen erally; 
the stellar mass derived assuming a iKroupal ([2001) IMF 
is only ~ 0.06 dex greater than M* from a Chabrier IMF, 
which is not able to account for the a values that we find 
(Table [lj. We note that based on our data alone we can- 
not distinguish whether this is due to a higher abundance 
of low- mass stars (Salpeter-like), or of higher-mass stars 
and their neutron star and black hole remnants. 

These results may appear to be inconsistent with those 
of iCappellari et akl (|2006| ). who find that an IMF with a 
higher normalization than Kroupa (e.g., Salpeter) leads 
to stellar masses that are sometimes greater than the 
dynamical masses. However, we note that the tension 
is entirely due to fast-rotating or lower-mass galaxies, 
neither of which are typical of the SLA CS lenses (e.g., 
IBarnabe et all 120091: lAuger et alJ l2"010ah : the SAURON 
data are consistent with our fin ding that massive ETG s 
do not have bottom-light IMFs (|Cappellari et al.ll2006l ). 

4.2. Adiabatic Contraction and a Varying IMF 

The data also are opposed to strong adiabatic contrac- 
tion, as modeled by B86. However, we cannot currently 
distinguish between a halo that has undergone mild AC 
and a NFW halo. This is partially due to imposing the 
form o f the M v ; r -M* relation; if we use the iDutton et alJ 
(2010) parameters for Equation 1 we find that for a fixed 
IMF the G04 model is strongly favored over the NFW 
model (and B86 is still disfavored), although we are again 
unable to distinguish between the two halos when the 
I MF is allowed to va ry with mass. 

ISchulz et all (|2009l ) similarly found that a Kroupa IMF 
and adiabatically contracted halo or Salpeter IMF and 
no AC adequately fit the data, but t hey do not consider 
a ma ss-dependent IMF. Likewise, Uiang fc Kochanekl 
(|2007l) find that a Salpeter-like IMF and AC model is pre- 
ferred (they do not directly compare their stellar masses 
with SPS model stellar masses, but their derived B-band 
M*/L = 7.2 M Q /L© is consistent with the M*/L we 
find for the SLAGS lenses assuming a Salpeter IMF, e.g., 
lAuger et al.ll2010al ). Note, however, that they do not test 
for a mass-dependent IMF (their M»/L is independent of 
mass). 

4.3. Beyond Structural Constraints 

The structural decomposition into bulge+halo mod- 
els is useful but may not be sufficient to fully under- 
stand how the IMF varies as a function of galaxy prop- 
erties. We conclude by considering how non-structural 
constraints might be used to further discriminate be- 
t ween CDM and a non-un iversal IMF. 

INapolitano et al.1 (|2010T ) use stellar dynamics and SPS 
models to determine the DM fraction for a sample of 
ETGs and compare these with the predicted DM frac- 
tions from a suite of bulge+halo toy models. They find 
that a Kroupa IMF with AC or Salpeter IMF with no AC 
best represent their data, although they do not directly 
fit the models and do not allow for a mass-dependent 
star-formation efficiency. Nevertheless, they find an in- 
triguing trend in which the IMF does not depend on 
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Fig. 1. — Bulge (dashed), halo (dot-dashed), and total mass (solid) profiles for each of the models investigated. The fiducial model, a 
NFW halo with Chabrier IMF, is shown in grey in each panel for comparison, and the vertical lines indicate the effective radius (dotted) 
and Einstein radius (dot-dashed); the radii probed by the stellar velocity dispersion and WL d ata are illustrated in the upper left panel. 
Note t hat the total mass profiles are approximately isothermal within the effective radius (e.g., Koopmans ct al. 2006, 2009; Auger ct al. 
2010a). The Chabrier IMF models require large halo masses in order to provide enough projected mass within the Einstein radius to fit 
the SL data. 



mass but may depend strongly on age in the sense that 
older galaxies have more Salpeter-like IMFs. As noted by 
INapolitano et all (|2010l ). this is the opposite trend than 
is expected from our relations. The cause of this discrep- 
ancy is unclear but may be related to how the samples are 
cho sen. The SLACS l enses are at higher redshifts than 
the INapolitano et al.l ([2010) galaxies and are selected to 
have spectra that match old stellar templates; they would 
therefore only span ages between approximately 5 and 10 
Gyr, for wh ich the age-IMF t rend i s essentially flat. 

Likewise, IGraves k. Faberl 1)20101 ) use trends between 
mass and metallicity ([Fe/H], [Mg/H], and [Mg/Fc]) to 
investigate the role of the IMF on the central dark matter 
fraction. They show that the a-enhancement of galax- 
ies (e.g., the [Mg/Fe] ratio) is mass-dependent and could 
therefore signal a more top-heavy IMF in more massive 
galaxies such that the relative number of core-collapse 
supernovae increases in more massive galaxies. How- 
ever, invoking more supernovae yields to super-abundant 
metallicities compared to the data, and another mecha- 
nism m ust therefore be posite d in order to remove these 
metals (|Graves fc Faberll2010D . 



5. CONCLUSIONS 

We have demonstrated that, given our standard as- 
sumptions, the IMF of massive ETGs cannot be "light" , 
even in the presence of significant adiabatic contraction. 
Furthermore, the IMF is only consistent with being uni- 
versal if the central dark matter profile is mass depen- 
dent, as in the AC models presented here. It is clear that 
if one wants to preserve a light and universal IMF one 
has to abandon standard assumptions about dark mat- 
ter halos, like the NFW profile. However, there are two 
important caveats. First, better constraints on the star- 
formation efficiency must be obtained from the data in 
order to draw definitive conclusions about the role of a 
mass-dependent IMF relative to CDM halo contraction. 
Second, although our conclusions are robust with respect 
to a wide variety of changes in our assumptions, we have 
not exhausted all possible families of theoretical models. 

We will address these caveats, quantify the relationship 
between the central dark matter profile and the IMF in 
more detail, and quantify the conseq uences of our various 
assumptions in a forthcoming paper ([Auger et al.ll2010bl 
in preparation). 
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TABLE 1 
Best Fitting Parameters and Goodness of Fits 



Halo 



IMF log[ M.,0 / M Q ] log[ M vir , / M ] 



x 2 /d.o.f. (InP) 



11.8 



11.6 



11.4- 



11.2 



11.0 



10.8 



NFW 


Chab 


11.38 ±0.09 


11.32 ±0.24 






2.68 


-19.16 




Salp 


11.29 ±0.04 


10.38 ±0.06 






1.18 


266.54 




Free 


11.44 ±0.06 


10.52 ±0.07 


0.08 ±0.04 


0.03 ±0.03 


1.04 


284.16 


G04 


Chab 


10.97 ±0.04 


10.11 ±0.11 






1.66 


189.76 




Salp 


11.36 ±0.05 


10.43 ±0.05 






1.12 


268.98 




Free 


11.26 ±0.06 


10.34 ±0.06 


0.03 ±0.05 


-0.08 ±0.02 


1.08 


281.58 


B86 


Chab 


10.94 ± 0.04 


10.07 ±0.10 






1.42 


228.43 




Salp 


11.47 ±0.06 


10.54 ±0.06 






1.26 


234.00 




Free 


11.25 ±0.06 


10.33 ±0.07 


0.01 ±0.05 


-0.11 ±0.02 


1.13 


269.05 



Note. — The 'Free' IMF is from Equation [5] for 77 = 0, log a = for a Salpeter IMF and log a R; — 0.25 
(log a ss —0.2) for a Chabricr (Kroupa) IMF. The % is given for a model with the mean parameters listed in the 
table, although the interpretation of Ax is difficult due to the large number of degrees of freedom; (lnP)is the 
mean of the natural log of the posterior probability. It is also difficult to directly interpret A(lnP), but values of 
several tens clearly indicate a strong preference for one model over another. 
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Fig. 2. — Best-fit M v i r -M* relation for each of the bulge+halo 
models considered. Note that the models with a Chabrier IMF 
(solid lines) require substantially larger halo masses at fixed stellar 
mass in order to fit the SL data. The grey lines indicate contours 
of constant star-formation efficiency. 
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tours), or B86 halo (dotted contours); the inner (outer) contour 
encloses 68% (95%) of the probability. The SPS masses were de- 
termined assuming the same IMF (Salpeter) for each galaxy and a 
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